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Abstract 

Our Universe is full of regions where extreme physical conditions are real- 
ized. Among the most intriguing cases are the so-called magnetars: neutron 
stars with very dense cores and super-strong magnetic fields. In this paper I 
review the current understanding of the physical properties of quark matter at 
r~| | ultra-high density in the presence of very large magnetic fields. I will discuss 

the main results on this topic, the main challenges that still remain, and how 
they could be related to the physics of magnetars. 
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in 

cn ■ 1 Introduction 

CN 

CN ' The realm of high density QCD pertains to situations of very high baryon density. 

At very high baryon density (and low temperatures) baryons get so squeezed that 
they start to overlap, thereby erasing any vestige of structure. Since in this case 
the quarks get very close to each other, the phenomenon of asymptotic freedom 
ensures their interactions to become weaker and weaker when the density becomes 
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5_i ■ higher and higher. At densities of the order of 10 times the nuclear density the 



quarks will be so weakly interacting that they can exist out of confinement. What is 
particularly interesting about cold and dense quark matter is that the fundamental 
QCD interaction is attractive in the color antitriplet channel. Once the quarks are 
deconfined and fill out the available quantum states up to the Fermi surface, this 
attractive interaction triggers the formation of diquark pairs at the Fermi surface, 
thus leading to the phenomenon of color superconductivity (CS) (for a historical 
account and detailed discussion of CS see [1]). 

In nature the combination of the high densities and relatively low temperatures 
required for color superconductivity can be found in the interior of neutron stars, 
which are the remnant of supernova explosions. At the same time, it is well-known [2] 
that strong magnetic fields, as large as B ~ 10 12 — 10 13 G, exist in the surface of regular 
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neutron stars, while in the case of magnetars they are in the range B ~ 10 14 — 10 15 
G, and perhaps as high as 10 16 G [3]. Moreover, the virial theorem pf] allows the field 
magnitude to reach values as large as 10 18 — 10 19 G. To produce reliably predictions 
of astrophysical signatures of color superconductivity, a better understanding of the 
role of the star's magnetic field in the color superconducting (CS) phase is essential. 

In recent years, several works |5]-[9] have been dedicated to elucidate the influence 
of a magnetic field in the ground state of CS matter. These investigations have 
revealed a richness of phases [9] with different symmetries and low energy properties. 

In order to grasp how a magnetic field can affect the color superconducting pairing, 
it is important to recall that in spin-zero color superconductivity, although the color 
condensate has non-zero electric charge, a linear combination of the photon and one 
of the gluons remains massless [TDl [EE], so the condensate is neutral with respect 
to the Abelian charge associated with the symmetry group of this long-range gauge 
field. This combination then behaves as the "in-medium" (also called "rotated") 
electromagnetic field in the color superconductor. Since this combination acquires 
no mass, there is no Meissner effect for the corresponding "rotated" magnetic field 
and consequently, a spin-zero color superconductor may be penetrated by a rotated 
magnetic field B. Moreover, it is worth to notice that despite all the superconducting 
pairs are neutral with respect to this long-range field, a subset of them is formed by 
quarks of opposite rotated charges Q. The interaction of the charged quarks with the 
magnetic field gives rise to a difference between the gaps getting contribution from 
pairs formed by oppositely charged quarks and those getting contribution only from 
pairs of neutral quarks. One consequence of such a difference is the change of the 
gap parameter symmetry [5] . If the field is strong enough, it actually strengthens the 
pairing of quarks of oppositely rotated charge [5]. One can intuitively understand 
this considering that the quarks with opposite charges Q and opposite spins, have 
parallel (rather than antiparallel) magnetic moments, so the field tends to keep the 
alignment of these magnetic moments, hence helping to stabilize the pairing of these 
quarks. 

Besides changing the symmetry of the gap and consequently the low-energy physics, 
a magnetic field can lead to other interesting behaviors too. It can produce oscilla- 
tions in the gaps and the magnetization |12j . the Hass-Van Alphen effect. Moreover, 
when the field strength is of the order of the Meissner mass of the rotated charged 
gluons, these modes become tachyonic [Z]-[8]. The solution to this instability is the 
formation of a vortex state of gluons which in turn boosts the magnetic field, creat- 
ing a peculiar paramagnetic state. This magnetic-field induced gluon vortex state is 
known as the Paramagnetic CFL (PCFL) phase. 

The paper is organized as follows. Section 2 and 3 outline the effects of a magnetic 
field in the gap magnitude and structure for three- and two-flavor spin zero color 
superconductors. In Section 4 I briefly mention several questions derived from the 
results here presented, as well as the main standing problems in the field of color 
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superconductivity and how taking into account the magnetic field (whether external 
or induced by the color superconductor) we could address some of them and have as 
a byproduct a potential solution for some of the puzzles in the physics of magnetars. 

2 Three Quark Flavors in a Magnetic Field 
2.1 MCFL Symmetry 

Let us consider a model of three quark flavors in the presence of a magnetic field 
at high baryon density. This system was investigated in [5] in the context of a NJL 
model based on the one-gluon exchange interaction of QCD. The first important 
thing to notice in this case is that a magnetic field affects the flavor symmetries of 
QCD, as different quark flavors have different electromagnetic charges. For three 
light quark flavors, only the subgroup of SU(3)l X SU(3)r that commutes with Q, 
the electromagnetic charge operator, is a symmetry of the theory. Based on the above 
considerations, and imposing that in the presence of an external magnetic field the 
condensate should retain the highest degree of symmetry, one can propose [5] the 
following ansatz for the gap structure in the presence of a magnetic field 
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The above gap is based on the quark representation ip T = (sj, s 2 , S3, d\, d?, d 3 , m, 112, U3). 
Notice that we ignored the symmetric gaps (see [6] for the complete set of symmetric 
and antisymmetric gaps in the gap structure in a magnetic field). For the purpose of 
symmetry considerations, they are no relevant, as they do not break any symmetry 
that is not already broken by the antisymmetric gaps. However, the reader should 
be aware that the symmetric gaps are nonzero. In general they are smaller than 
the antisymmetric gaps, because they originate from a color- repulsive, rather than 
attractive interaction. Nevertheless, we call attention that at magnetic fields of the 
order or larger than the baryon chemical potential scale, the symmetric gap that gets 
contributions from pairs of charged quarks can be of the order or larger than the 
antisymmetric gap that only gets contributions from pairs of neutral quarks [6]. 
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The order parameter ([I]) implies the following symmetry breaking pattern: 

SU(3) color x SU(2) L x SU(2) R x U(1) B x U^\l) A - SU(2) color+L+R . (2) 

The U^~\1)a symmetry is connected with the current which is an anomaly-free linear 
combination of s, d and u axial currents [13] . The locked SU(2) corresponds to the 
maximal unbroken symmetry, and as such it maximizes the condensation energy. 
Given that it commutes with the rotated electromagnetic group U(l) em , the rotated 
electromagnetism remains as a symmetry of the MCFL phase. 

The phase described by the order parameter (CEJ) is known in the literature as the 
Magnetic CFL (MCFL) phase. It locks SU(2) left and right flavor transformations 
with SU(3) color transformations, similar to the CFL phase, so it also breaks the chiral 
symmetry of the original theory. The low-energy physics of the MCFL phase differs 
from the CFL one, as it is characterized by five instead of nine Nambu-Goldstone 
(NG) bosons. One of the NG bosons is associated to the breaking of the baryon 
symmetry; three others are associated to the breaking of SU(2)a, and another one is 
associated to the breaking of U^~\1)a- The propagation of light in the MCFL phase 
is also different from the CFL all of the five NG bosons in the MCFL phase 

are Q-neutral. 

2.2 MCFL Effective Action 

The effective action of the three-flavor quark system in the presence of a magnetic field 
can be obtained using a Nambu- Jona-Lasinio (N JL) Lagrangian with the four-fermion 
interaction abstracted from one-gluon exchange [TO] . 

The mean-field effective action for such a theory can be written as 

I B ($,i>) = J {g^ojW^Jjo]" 1 ^^)^)^) +^(+)(x)[Gf +)0 ]-\x,y)^ i+) {y) 

x,y 

+</>(_) (x) [Gf.j ] _1 (x, y )V>(-) (y) + V»( )c 0*0 [G (o) ] _1 (x, y)ip {0 )c (y) 
+^ {+)c (x)[G W0 ]- 1 (x,y)ij {+)c (y) +^ { - )c (x)[G'[_ )0 )' 1 {x,y)il; H c(y)] 

+ 2$(o)c(x) A+ (x,y)i J (o)(y) + h.c] + -$ {+)c (x)A + (x,y)4> { _ } (y) 

+^ { _ )c (x)A + (x } y)ip w (y) + h.c.]} , (3) 

where the external magnetic field has been explicitly introduced through minimal 
coupling with the Q— charged fermions. The presence of the field is also taken into 
account in the diquark condensate A + = 75 A, whose color-flavor structure is given 
by Eq.p. 

In ()3]) symbols in parentheses indicate neutral (0), positive (+) or negative (— ) 
Q— charged quarks. Supra-indexes + or — in the propagators indicate, as it is cus- 
tomary, whether it is the inverse propagator of a field or conjugated field respectively. 
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Then, for example, [G^J -1 corresponds to the bare inverse propagator of positively 
charged conjugate fields, and so on. The explicit expressions of the inverse propaga- 
tors are 

[Gfo)o] -1 (*. V) = i^d, -m± ^°}6 4 (x - y) , (4) 
[Gt + y>]-\x,y) = [i 1 ^ ) -m±^°}8 4 (x-y) , (5) 
[Gf^-^y) = [iynjr' -m±n 1 °]5\x-y) , (6) 

with 

nf) = id, ± eA, . (7) 

Transforming the field-dependent quark propagators to momentum space can 
be performed with the use of the Ritus' method, originally developed for charged 
fermions [Hj and later extended to charged vector fields [15] . In Ritus' approach the 
diagonalization in momentum space of charged fermion Green's functions in the pres- 
ence of a background magnetic field is carried out using the eigenfunction matrices 
E p (x). These are the wave functions of the asymptotic states of charged fermions in 
a uniform magnetic field and play the role in the magnetized medium of the usual 
plane-wave (Fourier) functions e ipx at zero field. 

The transformation functions E^\x) for positively (+), and negatively (— ) charged 
fermion fields are obtained as the solutions of the field dependent eigenvalue equation 

{Tl^. 1 )E^\x)=E^\x){ 1 -p^) ) (8) 

with p^ given by 

P (±) = (po,0,±yf2\eB\k,p 3 ) , (9) 

and 

i?! ±) W = E^ ) («. (io) 

a 

with eigenfunctions 

E£\x) = Mn^e-^^^D^g^), (11) 
where D n Aq^) are the parabolic cylinder functions with argument Q(±) defined by 

Q{±) = j2\sW\(x 1 ±p 2 /eB) , (12) 
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and index n(±) given by 

eB 1 

n{±) = n {±) (k,a) = k±—~-a--, n (±) = 0, 1, 2, ... (13) 

2\eB\ * 

k = 0, 1, 2, 3, ... is the Landau level, and a is the spin projection that can take values 
±1 only. Notice that in the lowest Landau level, k = 0, only particles with one of the 
two spin projections, namely, a = 1 for positively charged particles, are allowed. The 
normalization constant A/^ (±) is 



N n(±) = (47r|e£|) WV)! • (14) 
The spin matrices 5(a) are defined as 

5(a) = diag(£ CT i, 5 a -i, 5 al , 5 a -x), a = ±1 , (15) 
and satisfy the following relations 

5 (±) f = 5 (±) , 5 (±) 5(±)=5 (±) , 5 (±) 5 (=f) = 0, (16) 

7 II 5(±) = 5(±)7 11 , 7 ± 5(±)=5(T)7 ± - (17) 

In Eq. ( TT7|) the notation 7" = (7°,7 3 ) and 'j ± = ('J 1 , 'J 2 ) was used. 
The functions E^ are complete 



k 

and ort honor mal, 



£ / dp dp 2 dp 3 E^(x)E p ±) (y) = (2n)^(x-y) 



I Ef(x)E^(x) = (27r) 4 A k 5 w 5(p -p' )5(p 2 -p' 2 )5(p 3 -p' 3 ) (19) 
with the (4 x 4) matrix A k given by 



(5(a = sgn[eB]) for k= 0, 

( 7 /or fc > 0. v y 

In Eqs. (H1D-([I9D we introduced the notation E { p ±] (x) = 7 (^ ±) (x)) t 7 . 
Under the E p (x) functions, positively (ip(+)), negatively (ip<-)) charged fields trans- 
form according to 

i'(±)( x ) = H dp dp 2 dp 3 E { p ±) (x)ip {±) (p) , (21) 
k J 
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VW X ) = E / dp dp 2 dp^ {±) {p)Ef\x) . (22) 
One can show that 

h,(U {+) , ± /^ ) - m]E^\x) = E^(x)[^(p^ ± /i^ ) - m] , (23) 

and 

[7 M (n ( _ )/i ± /^ ) - = E^(x)Mp^ ± /iV) - m] . (24) 

The conjugate fields transform according to, 

A+)c( x ) =E / dpodp 2 dp 3 E { - ) (x)ip (+) c(p), (25) 

^(_)c(x) =E/ dpodp 2 dp 3 E^\x)^ )c (p) . (26) 

After transforming to momentum space one can introduce Nambu-Gorkov fermion 
fields of different Q charges. They are the Q-neutral Gorkov field 



the positive 



and the negative one 



**> = ( t ) ' (27) 



*<+> ^ ( ) ■ (28) 



Using them, the Nambu-Gorkov effective action in the presence of a constant 
magnetic field B can be written as 

I B ®,1>) = \j ^%(p)S(o)(p)*(o)(p) 
+\ j -^^ {+) (p)S^(p)^ +) (p) + \ J , (30) 

where 

[Gt 0)0 ]-\P) A(o) \ 



5(0) (P) 



(31) 
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35 (p) 



5 ( -)(p) 



with 



'(+) L^(+)0J 

[G^ol-'b) A^_) 

A+ +) = fi_A+fi + , 
A+ } = fi + A+fi_, 

A(o) = ^oA + fioi 
Here we introduced the rotated charge-projector operators 

ft o = diag(l,l,0,l,l,0,0,0,l) , 
Q+ = diag(0, 0,0, 0,0, 0,1,1,0) , 
= diag(0,0, 1,0,0, 1,0,0,0) , 

which obey the algebra 

n v n v > = s w n v , rj, r{ = o, +, - . 

fi + ^+ + fi- = 1 • 

With the help of the charge projectors one can write the rotated char^ 
and the (0)-, (+/— )-charged fields as 

Q = ^ r]Vt ri = Q + — Q_ . 

r?=0,± 

and 

respectively. 



Notice that to form the positive (negative) Nambu-Gorkov field we used the posi- 
tive (negative) fermion field and the charge conjugate of the negative (positive) field. 
This is done so that the rotated charge of the up and down components in a given 
Nambu-Gorkov field be the same. This way to form the Nambu-Gorkov fields is 
mandated by what kind of field enters in a given condensate term, which in turn is 
related to the neutrality of the fermion condensate with respect to the rotated 

Q-charge. 

In momentum space the bare inverse propagator for the neutral field is 

[ G fo)o] _1 (p) = llM ± /"<W - m] , (44) 

where the momentum is the usual p = (po,Pi,P2,P3) of the case with no background 
field. 

For positively and negatively charged fields the bare inverse propagators are 

[Gf^ip) = MP^ ± /iM - m] , (45) 

and 

[G^ )0 ]-\p) = [iM^ ± ~ m] (46) 

respectively. 

2.3 MCFL at Strong Magnetic Fields 

The details of the gap equations obtained from ([3]) can be found in [6j. In the strong 
field limit eB ~ fi 2 , they become 



a r 9 2 f d 3 q A^ g 2 eB f A dq A^ 



3A 2 7a (2tt) 3 v /( g - /U )2 + 2 (A^) 2 3A 2 J -a (2tt) 2 ^ q _ ^ji + (A B) 

^ y rf 3 g /17 A A 7 A A x 

A ^4A 2 A(27r)3V 9 ^_^ )2 _ A 2 9^_^2 + 2 (Af)^ ' 1 




A* ~ 2 JA/z - /x 2 exp - , w -a\ • (49) 



with solution 



for the gap receiving contribution from pairs of charged quarks. It is instructive to 
look at the form of this gap. Just as in the conventional BCS gap, this strong-field 
MCFL gap goes as exp(— 1/NGj , where G = g 2 /3A 2 is the dimensionful coupling 
constant, and N represents the total density of states at the Fermi surface of those 
quarks contributing to the gap. In the absence of a magnetic field, the density of states 
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for a single quark is = [i 2 /2Ti 2 . The total density of states is then N = 4A^, as 
there are four quarks lending to each gap. When a magnetic field is present, the 
contributing quarks are shared among two terms, and Ng, depending on whether 
they are neutral or charged. Here, Ng = eB/4ir 2 . In the present case, we have 
iV = 2N fl + 2Ng, as two of the four quarks are charged. We then have the above 
expression fj49|) . 

On the other hand, the gap that only gets contribution from pairs of neutral 
quarks has solution 

A7 L / 36 21 1 3 v 

where x = g 2 fj 1 2 /A 2 ir 2 , and y = eB/fj, 2 and 

2^-2 \ 



AS" ~2jl^? exp [~^y (51) 

From ( H9l) it is clear that in the strong field limit, this gap increases with larger 
magnetic fields. If the strength of the field is such that eB > 2/i 2 , we see from 
comparing (jSTl) and (J151) that the MCFL gap surpasses the CFL gap, A^ > A^ FL - 



3 Two Quark Flavors in a Magnetic Field 

The case of two quark flavors in a magnetic field was recently considered in |16j . 
Since the quarks participating in the pairing are all charged with respect to the 
rotated electromagnetism, the external magnetic field does not change the structure 
of the gap, but only its magnitude. For 2SC pairing at zero temperature, the effective 
thermodynamic potential with an arbitrary magnetic field is given by 



O - O _i_ ^ ^ db 



eB 2eB 

E(2-M 



471-2 



n=0 



fjyJ^ e -2eBn-eBn In 



/i 2 — 2eBn + \i e 



V2eB 



n 



eB l2eB 



4^2 



E ( 2 - *no) 
n=0 



Pub 



fi 2 b — 2eBn — eBn In 



p 2 ub - 2eBn + [M ub 



V2eB 



n 



eB 00 ( 

" n=0 V 



8n0 



(p 2 +eBn\ 




p 2 + eBn + jj, ] + A 2 dp. 
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eB 00 / 



2tt 



n=0 



eB 00 ( 

Eh 



2tt 



n=0 



2 / 

2 j 



e 



e A2 



p| + e5n — ~p ) + A 2 + fyt 



p 2 + e5n — /i I + A 2 — 5^ 



dp 3 



^3. (52) 



where in order to impose the neutrality conditions and to satisfy (3- equilibrium 
constraints, we have to introduce all the chemical potentials for the conserved and 
commuting charges. The diagonal matrix of chemical potentials is given by 



Hij,a/3 — [H^ij — HeQij)^a/3 + ^/g^ 8 ^' (^s)a/9- 



(53) 



Here, /1, /x e , and /is are the quark, electron, and color chemical potentials respectively. 
The generators, Q and T 8 , are those of the electromagnetic group U(l) em and the color 
subgroup U(l)s- The individual quark chemical potentials then read 



Hur Hug 
Hdr Hdg 
Hub 
Hdb 



2 1 

H — + 3/^8) 

/i + g/i e + g//g, 

2 2 
A 1 — 3 He — 3/^8; 

1 2 
A 4 + 3/^ — 3/^8- 



(54) 
(55) 
(56) 
(57) 



By requiring 2SC quark matter to be invariant under the SU(2) C color subgroup, we 
avoid the need to introduce a second color chemical potential, 7/3, though in general 
this would be present, as there are two different color charges [17]. Following the 
notation of Ref . [T7] , we introduce the shorthand 



_ _ Hur + Hdg _ Hug + Hdr _ He H8 /ro\ 

h = g " 2 _// "T + T' (5bj 

r Hdg Hut Hdr Hug He (KC\\ 

OH = ~ 2 = ~ 2 = y • {^) 

For large magnetic fields, eB ~ 4/i 2 and assuming the density to be large enough 
to avoid gapless modes, the thermodynamic potential ( 152]) reduces to 



n 



A 2 



h%> 

12n 2 



eB 
I^ 2 



HI 



eB 



Hub 



eB 
2^ 



(vW7l) 2 + A 2 + ^( P3 -Jlf + A^ 



2 J 
dp 3 . 



(60) 
(61) 
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The 2SC gap that solves the gap equation derived from this strong field potential 



is 



A = 2^A 2 - 



/x 2 exp 



■—) 

2GeB 



(62) 



In 2SC, all four of the participating quarks carry a charge. The total density of 
states is then 4[(e/2)S/47r 2 ] = eB/2n 2 . Defining G = 4G, the exponential becomes 



A ~ exp (-gTb) 



(63) 



which is the usual BCS form. 

Imposing color and electric neutralities one can show that the solutions for the 
chemical potentials in the strong field limit are 



3 5 R 



+- 



2 (-12eBfi + VQ^/(eB) 3 + 24(e£) V 



3(-12e£// + VQj(eB) 3 + 24(eB) 2 fi 2 



3/x - 2/x 8 



(64) 
(65) 



4 Final Remarks 

This paper has mainly outlined the effects of a magnetic field in the gap magnitude 
and structure of three and two quark flavor spin zero color superconductors. Some 
open questions naturally emerge from these results. First, it will be interesting to 
investigate how the smaller number of NG bosons in the MCFL phase and the fact that 
none of them are charged is reflected in the comparison of the transport properties 
of MCFL and CFL. How will the gap parameter and chemical potentials behave at 
lower magnetic fields in the two-flavor case? How do the free energies of the two 
and three flavor systems compare when in addition to the magnetic field a small, but 
finite temperature is also introduced? 

However, the most urgent question in the field of color superconductivity is related 
to the stable ground state that realizes at intermediate densities. This problem arises 
at intermediate densities when the Fermi surface imbalance between pairing quarks 
becomes too large thus giving rise to gapless modes. The imbalance itself comes 
from imposing beta equilibrium and neutrality or when the density is such that the 
strange quark mass cannot be ignored any longer. As it turns out, once the gapless 
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modes occur, there are also chromomagnetic instabilities present, indicating one is not 
working in the correct ground state. A solution to this chromomagnetic instability 
on which an inhomogeneous gluon condensate and an induced magnetic field are 
generated hence modifying the ground state of the superconductor was proposed in 
[IB] . Since it involves the spontaneous generation of a rotated magnetic field by the 
gluon condensate, it could also serve to address some of the puzzles of the standard 
magnetar model [19], thereby connecting magnetars and color superconductivity. 

On the other hand, a magnetic field can have some other effects on the color 
superconductor that have not been explored yet and which could give rise to unex- 
pected and interesting new physics. One possibility is that the Dirac structure of 
the gap may also include an anomalous magnetic moment term that is formed by 
the contraction between the spin operator and the field. If such a term breaks the 
same symmetry as the gap with Dirac structure C75, it has all the right to be present 
once the gap is generated, since it is not protected by any symmetry after that. This 
reasoning is similar to the impossibility to avoid the presence of symmetric gaps in 
spin zero superconductivity. The need to consider a dynamical anomalous magnetic 
moment along with a dynamical mass was recently studied in the context of massless 
QED and magnetic catalysis of chiral symmetry breaking [20]. Over there it gave 
rise to a nonperturbative Zeeman effect. It would be worthy to find out what conse- 
quences a nonperturbative anomalous magnetic moment term could have on the color 
superconductor and particularly, whether it can offer a different kind of solution to 
the intermediate density instability problem. 

I am grateful to the organizers of the CSQCD conference for their invitation. 
This work was supported in part by the Department of Energy Nuclear Theory grant 
DE-FG02ER41599. 
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